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Abstract 

r — > Trvrvrz^T- decays are analysed within the framework of the resonance effective 
theory of QCD. We work out the relevant Lagrangian that describes the axial-vector 
current hadronization contributing to these processes, in particular the local ai(1260)- 
p(770)-Goldstone interactions. The new coupling constants are constrained by impos- 
ing the asymptotic behaviour of the corresponding spectral function within QCD. Hence 
we compare the theoretical framework with the experimental data, obtaining a good 
quality fit from the ALEPH spectral function and branching ratio. We also get values 
for the mass and on-shell width of the ai(1260) resonance. In this way we are able to 
provide the structure functions that have been measured by OPAL and CLEO-II and 
we find an excellent agreement. 
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1 Introduction 



The decays of the r lepton offer an excellent laboratory for the analysis of various topics in 
particle physics. In particular, r decays into hadrons allow to study the properties of the 
vector and axial-vector QCD currents, and provide relevant information on the dynamics 
of the resonances entering into the processes. The hadronization of those currents involves 
the strong interaction at low energies and, therefore, non-perturbative features of QCD have 
to be implemented properly into their evaluation. The Effective Theory framework has a 
long-history of successful achievements in this task. 

At very low energies, typically E <C Mp (being Mp the mass of the p(770), the lightest 
hadron resonance), chiral perturbation theory (xPT) 0|2] is the corresponding effective 
theory of QCD. However the decays r ttttttz/t-, through their full energy spectrum, happen 
to be driven by the p(770) and ai(1260) resonances, mainly, in an energy region where they 
reach their on-shell peaks. In consequence xPT is not directly applicable to the study of 
the whole spectrum but only to the very low energy domain j3|. Until now the standard 
way of dealing with these decays has been to use O(p^) xPT to fix the normalization of 
the amplitudes in the low energy region and, accordingly, to include the effects of vector 
and axial-vector meson resonances by modulating the amplitudes with ad hoc Breit-Wigner 
functions 0113] • However we will see that this modelization is not consistent with O(p^) xPT, 
a fact that could potentially spoil any outcome provided by this procedure. 

In the last years several experiments have collected good quality data on r ^ vrTrvrz/^, 
such as branching ratios and spectra jHllZj or structure functions jB] , and their analysis within 
a model-independent framework is highly desirable if one wishes to collect information on 
the hadronization of the QCD currents. The Effective Field Theory still provides the proper 
tool to proceed in this study. At energies E ~ Mp the resonance mesons are active de- 
grees of freedom that cannot be integrated out, as in xPT, and they have to be properly 
included into the relevant Lagrangian [H]. The procedure is ruled by the chiral symmetry 
under S'[/(3)l (S> SU{3)b., that drives the interaction of Goldstone bosons (the lightest octet 
of pseudoscalar mesons), and the SU{3)y assignments of the resonance multiplets. Its sys- 
tematic arrangement has been put forward in Refs. |10| [TT] as the Resonance Chiral Theory 
(RxT). We will attach to this framework and we will complete it, when needed, in order to 
fulfill our requirements. 

A complementary tool, thoroughly intertwined with the R^T, is the large number of 
colours (Nc) limit of QCD. Strong interactions in the resonance region lack an expansion 
parameter that, as in xPT, could provide a perturbative treatment of the amplitudes. How- 
ever it has been pointed out ^2] that the inverse of the number of colours of the gauge group 
SU{Nc) could accomplish this task. Indeed large- Ai'c QCD shows features that resemble, 
both qualitatively and quantitatively, the Nc = 3 case [12]. Relevant consequences of this 
approach are that meson dynamics in the large-A^^c limit is described by the tree diagrams 
of an effective local Lagrangian; moreover, at the leading order, one has to include the con- 
tributions of the infinite number of zero-width resonances that constitute the spectrum of 
the theory. 



1 



In this article we study the r — > innri'T processes within the Resonance Effective Theory 
of QCD. We will evaluate the form factor associated to the — 1+ piece of the axial-vector 
current that is, by far, the dominant contribution to these processes. To proceed we will work 
out the effective Lagrangian that endows the dynamics of the decay at the leading order in 
the 1/Nc expansion, though we will consider the lightest vector and axial- vector octets of 
resonances only. Moreover we will provide these with a finite width, due to the fact that 
they do really resonate along the full energy spectrum. In order to improve the construction 
of the Effective Theory we will constrain the new unknown couplings in the Lagrangian by 
demanding the asymptotic behaviour of the form factor ruled by perturbative QCD. The 
theoretical construction will be followed by the phenomenological analysis of experimental 
data. 

The article is organized as follows : The resonance effective theory is studied in Sect. 
2. In Sect. 3 we consider the r — > ttttt^Vt processes and calculate, within our approach, the 
relevant hadronic matrix elements. Sect. 4 is devoted to discuss the constraints imposed 
by perturbative QCD on the form factors driven by the resonances, which lead to a set 
of relations between the coupling constants. The analysis of experimental data within our 
framework is presented in Sect. 5, while in Sect. 6 we sketch our conclusions. An Appendix 
recalls the basic features of the very-low energy domain of the t~ — > 7r'^7r'"7r~i/,- process 
within xPT. 



2 The Resonance Effective Theory of QCD 

The very low-energy strong interaction in the light quark sector is known to be ruled by the 
SU SU {?>)y{, chiral symmetry of massless QCD implemented in xPT. The leading O(p^) 
Lagrangian is : 

— {uX + X+), (1) 



r(2) 



where 



X± = u^XU^ ± ux^u , X = 2So(s + ip) , (2) 
and (...) is short for a trace in the flavour space. The Goldstone octet of pseudoscalar fields 
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a=l 



V 



TT 

K' 



7r+ 



V2 



(3) 



is realized non-linearly into the unitary matrix in the flavour space 

^ 1 



u{i~p) — exp < i 



V2F 



(4) 
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that transforms as 

9r u{ip) h{g, ip^ = h{g, ip) u{ip) g{ , (5) 

with g = {gL,gR) £ SU{3)i, ® SU{3)b. and h{g,ip) G SU{3)v- External hermitian matrix 
fields r^, i^, s and p promote the global SU (3)l ® SU (3)r symmetry to a local one. Interac- 
tions with electroweak bosons can be accommodated through the vector = (r^ + £^)/2 and 
axial-vector = {r^ — £^)/2 fields. The scalar field s incorporates explicit chiral symmetry 
breaking through the quark masses s = M. + . . ., with M. = diag(m„, m^, m^) and, finally, 
F ~ ~ 92.4 MeV is the pion decay constant and BqF^ = —{0\iIjiIj\0)o in the chiral limit. 

The final hadron system in the r ttttttUt- decays spans a wide energy region 3m^ ^ 
E < Mr that is populated by resonances. As a consequence an effective theory description 
of the full energy spectrum requires to include the resonances as active degrees of freedom 
and, therefore, to consider R^T. This recalls the ideas of effective Lagrangians from Ref. 0. 
The interaction of the Goldstone bosons (the lightest octet of pseudoscalar mesons) with the 
resonances is driven by chiral symmetry, while the resonances enter as SU (3)v multiplets. In 
Ref. jini the simplest ^ couplings were introduced using the antisymmetric tensor formulation 
to describe the vector and axial-vector resonances but keeping linear terms in the latter, i.e. 
allowing for the couplings of the pseudoscalars to only one resonance. This formalism proved 
successful in accounting, upon integration of resonances, for the bulk of the low-energy cou- 
pling constants at in xPT with no need of the addition of local terms ^T]. We will 
attach to this antisymmetric tensor formulation and, accordingly, we will not include the 
0{p'^) xPT couplings into our effective action in order to avoid double counting of contribu- 
tions. The success of RxT in the description of many observables has been remarkable, in 
particular for the pion electromagnetic form factor entering r~ — >■ tt'tt^Uj. 

Therefore in order to introduce the octet of resonance fields we define III 



R^^E>^aRa, (6) 
V2 a=l 

where Ra = Va, Aa, stands for vector and axial-vector meson fields, respectively, that trans- 
form as 

R -> hig,^)Rh{g,^y (7) 

under the chiral group. The fiavour structure of the resonances is analogous to the Goldstone 
bosons in Eq. (jH)). We also introduce the covariant derivative 

V^X = d^X+[T^,X] , (8) 
= -[u'^{df,-irf,)u + u{d^-i£i,)u'^] , 

for any object X that transforms as R in Eq. ((Tj), like and x±- 

^ These were including one 0(p^) chiral tensor in the interaction operators. It is not possible to construct 
an interaction term with lower order in momenta. 
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Let us now restrict the analysis to the case of V{1 ) and y4(l++) resonances, in order to 
describe the couphngs of p(770) and ai(1260) states, respectively. Moreover we include the 
lightest SU (3)v octet of both vector and axial-vector resonances only. Although the leading 
1/Nc expansion requires to consider an infinite number of resonances their contribution 
happens to be suppressed as their masses become heavier. 

The kinetic terms for the spin 1 resonances in the Lagrangian read, in the antisymmetric 
tensor formulation, 

1 

= --{V^Rx.V^R'^ - -fR,uRn ^ R = V,A, (9) 

being My, Ma the mass of the octet of vector and axial-vector resonances, respectively. The 
lowest order interaction Lagrangian includes three coupling constants |10j. 



= ^(V/-"'^), (10) 

where f^" = uFi^u^ ± F^u and -Fr,l are the field strength tensors associated with the 
external right and left fields. All coupling parameters Fy, Gy and Fa are real. Notice 
that the 0{p^) interaction Lagrangian (fTUI) does not allow a tree- level coupling of the form 
ai(1260)7r7r7r. 

In the case of r ^ tttttti't-, in addition to one-resonance-exchange diagrams one has to 
consider the contribution given by the chain ai(1260) — >■ p(770)7r — tttttt, which is mediated 
by both vector and axial-vector resonances. Thus we need to go one step beyond the analysis 
in Ref. including bilinear terms in the resonance fields that lead to a coupling of the 
form aipn. For the processes under study, we only need to consider terms that include a 
vertex with one pseudoscalar. This goal is achieved with one 0{p^) chiral tensor. Hence the 
most general interaction Lagrangian allowed by the symmetry constraints can be written as 

£rP = EA.OVAP, (11) 

i=l 

where A, are new unknown real adimensional couplings, and the operators CyAP given by 

(12) 



Wap 


= ( 




) , 


'-^VAP 


= i 




ax 

fi 1 ' 


*-^VAP 


= i 


([V%,, A'^" 




'-^VAP 


= i 






^VAP 


= i 


([V"V, A'^'^ 


] ) 



with h^y = V fj,Uy + V yU^. We emphasize that this set is a complete basis for constructing 
vertices with only one pseudoscalar; for a larger number of pseudoscalars additional operators 
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may emerge. As we are only interested in tree level diagrams, the equation of motion arising 
from ;^PT, 

K = ^x-- ^(X-), (13) 

has been used in C^^^ in order to eliminate one of the possible operators. 

In summary we will proceed in the following by considering the relevant Resonance Chiral 
Theory given by : 

/:kxt = + cl^ + ct^ + ^ + c^ + £r • (14) 

It is important to notice that >Cr^t is not an effective theory of QCD for arbitrary couplings 
in the interaction terms C^^^ and C^^^ . These carry some symmetry properties of QCD, 
implemented in the operators, that do not give any information on the coupling constants. 
However we will see in Section 4 that well accepted dynamical properties of the underlying 
theory provide some information on those couplings. This is as far as we can go without 
including modelizations. 



3 Axial— vector current form factors in r ^ ttttttz/ 



r 



In this Section we describe the r — > ttttttz/t- decays in the isospin limit and proceed to calculate 
the relevant hadronic matrix element within the resonance effective theory given by Z^r^^t- 
The decay amplitude for the n^n^n^Ur and — ^ n^n^n^Ur processes can be written 

as ^ 

M± = -^KdM..7^(l-75)wrT±^ , (15) 

where V^d is an element of the Kobayashi-Maskawa matrix and T-j.^ is the hadronic matrix 
element of the participating — QCD currents. In the isospin limit there is no contri- 
bution of the vector current to these processes and, therefore, only the axial-vector current 
appears : 

T±,.(Pi,P2,P3) = (7ri(pi)7r2(p2)7r^(P3)| AXO)|0) , (16) 

being vr^ the one in r~ — h^tt'u'Vt- and vr" that in r~ — >■ ir'n^TT^Vr (in the following we 
will omit the subscript ± when speaking generally). The hadronic tensor can be written in 
terms of three form factors, Fi, F2 and Fp, as [1^] 

= Vl'F^ + V^F2 + Q^Fp , (17) 



where 



= p^+p^+p^ . (18) 
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Now if r stands for the r 
written as 

|2 



-r-r-rVt- decay width, the spectral function dT/dQ^ can be 



where s ■ 
given by 



dT 

= {Pi 



ud 



128 {2n)^Mr 




ds dt 



(19) 



P3) 



t = {p2 +^3)^5 and the hadronic structure functions Wsa and Wa are 



WsA 

Wa 



iQ'' Fp) {Q, F*] 



(20) 



The phase-space integrals extend over the region allowed for a three-pion state with a center 
of-mass energy \/Q^: 

I ds 

lAml 



dsdt = 



t+{s) 



dt 



where 



,2\2 



(21) 
(22) 



1 

with A(a, 6, c) = (a + 6 — c)^ — Aab. We have neglected here the mass of the Vr- 

In the decomposition of in Eq. (fTTjl . the form factors Fi and F2 have a transverse 
structure in the total hadron momenta and drive a = 1"^ transition. Bose symmetry 
under interchange of the two identical pions in the final state demands that Fi{Q'^,s,t) = 
F2{Q'^,t,s). Meanwhile Fp accounts for a. = 0^ transition that carries pseudoscalar 
degrees of freedom. As a consequence, the conservation of the axial-vector current in the 



chiral limit imposes that Fp\ 



0, and the scalar form factor must vanish with the square 



of the pion mass. Hence its contribution to the decay processes will be very much suppressed. 



3.1 Evaluation of the matrix amplitude 



We proceed now to calculate the hadronic amplitudes as given by our Resonance Effective 
Theory described in Section 2. Within the large- A^^c framework, one should evaluate all tree- 
level diagrams generated by -Cr^t that contribute to the decays ^. 

In the low region, the matrix element in Eq. p6|) can be calculated using xPT. At 
0{p'^) one has two contributions, arising from the diagrams (a) and (b) of Fig. ^ The sum 
of both graphs yields 



± 



2^2 

2V2 f 



9 HP 



QfiQu 



(2p3 -p2 -PiT - 



K± ml 



mt 



— S 
3F \ 



2fi 



ml [3 {u-ml)~Q\l±2K^)] 

2Q2(Q2_^2) 



Q 



(23) 



^We remind that though large- A^c enforces the contributions of the infinite spectrum of QCD resonances, 
we only include the lightest octet of vector and axial-vector mesons. The phenomenology seems to support 
well this approach. 
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(b) 



Figure 1: Diagrams contributing to the hadronic amplitude T±^ at in xPT. The crossed 

circle indicates the insertion of the axial-vector current. 



where we have defined u = {pi + = Q"^ — s — t + 3m^, k+ = 1, = 1/2. As expected 
from PCAC, it can be seen that the amplitude is transverse in the limit where the pion mass 
is neglected. 

In the limit of low (i.e. <^ M|.) the resonance fields can be integrated out, and 
their contributions should reduce to the corrections to Eq. (I23|) obtained from the effective 
0{p'^) interactions and higher order couplings in the standard xPT Lagrangian. In hadronic 
T decays, however, we need to implement the resonance degrees of freedom explicitly because 

~ M|.. As stated in Section 2 we assume here that the resonances p(770) and ai(1260) 
give the dominant contributions to r — vrvrvrz/^ decays. Consequently, in addition to the 
usual xPT diagrams leading to Eq. we include resonance-mediated contributions to the 
amplitude, to be evaluated through the interacting terms C^^^ and C^^^^ of the resonance 
effective theory. 

The relevant diagrams to be taken into account are those shown in Fig. |21 Let us start by 
considering the contributions given by the graphs in Figs. El (a) and (b), which involve only 
one intermediate resonance. From the interaction Lagrangian in Eq. ()10|). the sum of both 
graphs leads to 



(24) 



3F3 

where we have defined 



ai{Q^, s,t) 



+ 



,s,t) 



s 

3 



V 



j^2Gy 



2Q^ - 2s 



u 



V 



M, 



+ 



u 



t-Mt 



u — t 



mt s 



Mt 



(25) 



Fy Q Q — IIU^ S — IVly 

As expected, the pseudoscalar form factor a2(<5^, s,t) is found to vanish in the chiral limit. 

The two-resonance contribution in Fig. |2l (c) can be evaluated taking into account the 
chiral couplings from C^^^ in Eq. ()12|1 . We obtain : 



T. 



{2R) 



AFaG 



V 



3 F3 - MX 



(3{Q\s,t)V,^ + (3{Q\t,s)V, 



2n 



(26) 
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(a) (b) (c) 

Figure 2: Resonance-mediated diagrams contributing to the hadronic amplitude T^^ , (a) 
and (b), and T^^\ (c). 

where 

m\ s,t)^~i (v + A") ^ + F(Q\ s) ^±±^ + t) ^ . (27) 



with 



Ao 



1 

72 



4Ai +A2 + Y + A5 



A' = 71 (^2-A3 + ^ + A5. 



A" = ;^(A2-y-A5). (28) 

Notice that the form factor f3{Q'^,s,t) is written in terms of only three combinations of the 
unknown couplings in C^^^, namely Aq, A' and A". 

The final result from our evaluation gives, accordingly, the addition of all the amplitudes 
considered here : 

= + Ti;^) + Tif^ , (29) 

and, for later use, if we only consider the = 1+ mediated form factors Fi and F2 we will 
specify T±^ 

3.2 Implementation of off- shell widths 

The form factors in Eqs. (^3)) and ^T7\ include zero-width p(770) and ai(1260) propagator 
poles, leading to divergent phase-space integrals in the calculation of the r —>■ TrnnUr decay 
width as the kinematical variables go along the full energy spectrum. The result can be 
regularized through the inclusion of resonance widths, which means to go beyond the leading 
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order in the I/Nq expansion, and implies the introduction of some additional theoretical in- 
puts. Moreover both contributing resonances, especially the ai(1260), are rather wide. Hence 
energy-dependent widths should be included in order to handle their off-shell character. This 
issue has been analysed in detail within the resonance chiral effective theory in Ref. JHl- For 
the case of spin-1 vector resonances, it is seen that one can define the off-shell width by 
taking into account the pole of the two-point function of the corresponding vector current 
that arises from the resummation of those diagrams that include an absorptive contribution 
of two pseudoscalars in the s channel. The width is defined then as the imaginary part of 
this pole. In the case of the p(770) meson this procedure leads to the result 



Mv Tp{s 



96 ttF^ 



1 



(30) 



where ap = yl — 4mp/s, and 9{x) is the step function. In the case of the ai(1260) meson 
the construction following the definition above, though well defined, is much more involved. 
It would amount to evaluate the axial-vector-axial-vector current correlator with absorptive 
cuts of three pions within the resonance effective theory. This corresponds to a non-trivial 
two-loop calculation within a theory whose regularization is still not well defined. This 
calculation is far beyond our scope. In order to perform the phenomenological analysis we 
will introduce a chiral-based off-shell width for the ai(1260) resonance that endows the 
appropriate kinematical and dynamical features : 

r., (Q^) = r., (Ml) ^ (^^J 9{Q'-9 ml) , (31) 



where 



J dsdt i^V^ \BWp{s)\^ + Vi \BWp{t)\^ + 2 (Fi ■ V2) Re[BWp{s) BWpit)*] 



(32) 



and 



BWpis) = ' (33) 

My — s — z MvTp{s) 

is the usual Breit-Wigner function for the p(770) meson resonance shape, and the energy- 
dependent width Tp{s) is given by Eq. (jHUIl . 

In conclusion, while the off-shell width of the p(770) meson does not introduce any ad- 
ditional free parameters on our result for T±p, the ai(1260) width includes two new ones, 
namely the on-shell width Ta^lM^) and the parameter a ruling the behaviour in Eq. ()31|). 

4 Asymptotic behaviour and QCD constraints 

As commented at the end of Section 2, our Lagrangian Cr^t^ -^q. (fT^ does not provide 
an effective theory of QCD for arbitrary values of its couplings. The only tool used in its 
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construction were relevant symmetry properties of QCD, but it is clear that the full theory 
should unambiguously predict the values of the coupling constants. Although it is not known 
how to achieve this goal from first principles, several ideas based on matching procedures 
have been developed ^TJEl- In this Section we show how to proceed in order to obtain 
information on the, a priori, unknown couplings in C^^^ . 

The QCD ruled short-distance behaviour of the vector form factor in the large- A'c limit 
(approximated with only one octet of vector resonances) constrains the couplings of £3 
Eq. (fTUj). which must satisfy |TTj : 

1 - ^ = . (34) 



In addition, the first Weinberg sum rule |18j, in the limit where only the lowest narrow 
resonances contribute to the vector and axial-vector spectral functions, leads to 

F^-Fl = F\ (35) 

In Ref. ^] it was also noticed that there is an additional constraint coming from the axial 
form factor, namely : 

2FvGv - = , (36) 

that is very well satisfied by the phenomenology. However this relation is not necessarily 
true when we consider the inclusion of vector-pseudoscalar-axial couplings as given by C^^^ . 
In the following we are going to enforce the relation given in Eq. (jHUj) and we will see what 
happens when we relax this last condition. From Eqs. ()34l35l36p . the three couplings Fy, 



Gv and Fa in Eq. (fTUj) can be written in terms of the pion decay constant : Fy = a/2F, 
Gv = F/ \pl and Fa = F. We will adopt these results throughout this paper, unless stated 
otherwise. 

An analogous analysis should be done in the case of the axial two-point function Ua'IQ^), 
which plays in r — Trvrvrz/,- processes the same role than the vector-vector current correlator 
does in the r n-RVr decays, driven by the vector form factor. The goal will be to obtain 
QCD-ruled constraints on the new couplings of C"^^^ , Eq. (fTTjl . similar to those obtained 
above. As these couplings do not depend on the Goldstone masses we will work in the chiral 
limit but our results will apply for non-zero Goldstone masses too. In the chiral limit the 
nT(Q^) correlator becomes transverse, hence we can write 

nr(Q') = i.Q'Q'' - d^^'Q^) nA(Q') • (37) 

As in the case of the pion and axial form factors, the function n^((5^) is expected to satisfy 
an unsubtracted dispersion relation. This implies a constraint for the J = 1 spectral function 
Imn^((5^) in the asjTiiptotic region, namely 

Imn^(Q2) , . (38) 
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Now, taking into account that each intermediate state carrying the appropriate quantum 
numbers yields a positive contribution to Imn^(Q^), we have 



Imn^(Q2) > 



3Q2 



(39) 



(i$ being the differential phase space for the three-pion state. The constraint in Eq. 
then implies 

Wa 



lim 

32^00 Jo 



ds 



dt 



2^2 







(40) 



where Wa is the structure function defined in Eq. (j2(jp . It can be seen that the condition in 
Eq. (pn|) is not satisfied in general for arbitrary values of the coupling constants in the chiral 
interaction Lagrangian. Indeed, once Fy, Fa and Gy have been fixed by Eqs. (jHH), and 
()36|) . we find that the constants A' and A" appearing in the two-resonance contribution to 
the hadronic tensor are required to satisfy the relations 



2 A' 



- 1 
A" 








(41) 



Now we can come back to our results for T-j.^ in Section 3 and enforce the QCD driven 
constraints on the couplings of our Effective Theory as given by Eqs. ()34|35|36|41|) . We 
have : 



1+ 



^(IR) 



1+ 



2^2 



2s- Ml 



^1. + 1 - 



t 



2 t - Ml 



2n 



(42) 



and 



where 



/3(Q^s,t) 



± 



2V2 



Ml 



(3{Q\s,t)Vi^ + (3{Q\t,s)V, 



2fj. 



2 s-M^ 



2Q' 



Ao 



(43) 



(44) 



Thus, it turns out that the hadronic amplitude can be written in terms of only one unknown 
coupling parameter, namely Aq. 

As mentioned in Section 2, the resonance exchange approximately saturates the phe- 
nomenological values of the 0{p'^) couplings in the standard xPT Lagrangian. This allows 
to relate both schemes in the low energy region, and provides a check of our results in the 
limit <C My. We have performed this check (see Appendix A), verifying the agreement 
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between our expression Eq. (jl^ and the result obtained within xPT in Refs. jSlEj coming 
from saturation by vector meson resonances of the 0{p'^) couphngs : 



rpxPT 



2V2 



1 + 



3s 



3t 



2fj. 



+ chiral loops + 0{p^) . (45) 



As an aside, it is worth to point out that this low-energy behaviour is not fulfilled by all 
phenomenological models proposed in the literature. In particular, in the widely used model 
by Kiihn and Santamaria (KS) j4j the hadronic amplitude satisfies 



(KS) 



1 + 



t 



Ml, 



2^1 



(46) 



Thus, while the lowest order behaviour is correct (it was constructed to be so), it is seen that 
the KS model fails to reproduce the xPT result at the next-to-leading order. Accordingly 
this model is not consistent with the chiral symmetry of QCD. 



5 Phenomenology of r ^ nnnVr processes 

In this Section we discuss the ability of our Effective Theory to describe the experimental 
observations in r ^ ttttttz/^ decays, taking into account the data obtained from the measure- 
ments of r ^ -R-K-RVj. spectra, branching ratios and structure functions. 

Before proceeding let us specify which are the parameters left unknown in our evaluation 
of the hadronic matrix amplitude. To analyse the experimental data we will only consider the 
dominating = l"*" driven axial-vector form factors, being the pseudoscalars suppressed by 
pion mass factors. Hence we notice that ^+^11+ = — T_^|i+ and, accordingly, we have the 
same predictions for both tt+tt^tt^z/^ and -k^ti^tx^Vt- processes in the isospin limit. 

As we have shown in the previous Section, the requirement of the proper asymptotic behaviour 
of vector and axial-vector spectral functions consistent with QCD imposes several constraints 
on the couphng constants of our Lagrangian £r;^t in Eq- (fT^ and, as a consequence, the 
hadronic amplitude for the decays r -r-r-rVj. only involves one adimensional unknown 
combination of coupling constants, Aq. In addition, we have three remaining parameters, 
related with the ai(1260) resonance : these are its mass M4, its on-shell width Va^{M'j^ and 
the exponent a, which has been introduced in Eq. ()31|) to account for the dominating energy 
dependence of the off-shell ai(1260) width. The mass of the vector octet My, to be identified 
with the p(770) mass, is better obtained from the vector form factor of the pion ^JI^O] and 
we take here My = 775.1 MeV. In this way, we deal with four free parameters to fit the 
experimental data. 

5.1 Fit to the ALEPH branching ratio and spectral function 

In order to check the consistency of the available experimental results with the theoretical 
description provided by our chirally-based resonance framework we have carried out a fit of 
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the unknown parameters entering the amphtudes. We have chosen to fit the experimental 
values for the tt+tt^tt^z/^ branching ratio and normalized spectral function obtained by 

ALEPH jH], which quotes the unfolded spectrum for the decay including the corresponding 
covariance matrix. The data are collected in 57 equally spaced bins with a Svr squared 
invariant mass ranging from = 0.275 GeV^ to = 3.075 GeV^. We have noticed that 
bins n = 2 and n = 54 produce anomalously high contributions, due to their tiny errors, 
and we have discarded them. In addition we have fitted the total branching ratio to the value 
quoted by ALEPH [HI, BR(r^ ti^-k^-k^Ut) = (9.15 ± 0.15)%, which has been introduced 
as an additional fitting point. In order to control possible fake results the procedure has 
been carried out both with the MINUIT package j2I] and with an independent minimization 
procedure. As output of our fit we find the four parameter set 

Ao = 11.9 ±0.4 , 

a = 2.45 ±0.15 , 

Ma = (1.204 ± 0.007) GeV , 

Ta,{Ml) = (0.48 ±0.02) GeV , (47) 

with x^/d.o.f. = 64.5/52 (these numbers include only statistical errors). We find that the 
values of the parameters in Eq. (fTTjl lead to BR(r^ — n^n^n^Ur) = (9.1 ± 0.3)%. The 
theoretical spectrum arising from the central values of our fit, together with ALEPH data, 
are represented in Fig. El It can be seen that the agreement is good, if one considers the 
very small errors in the experimental data. On the other hand, the values obtained for the 
ai(1260) mass and width are in good agreement with the world average values quoted by the 
Particle Data Group 22 . 

Taking into account the value of a in Eq. (j47|) . we have carried out three additional three- 
parameter fits in which a has been fixed to the values 2, 2.5 and 3, respectively. The results 
are shown in Table I. Although the results clearly prefer a value of a = 2.5, one can see that 
the quality of the fits is reasonably similar to that of the four-parameter case, and the values 
of Ao, Ma and {M'^) are not significantly modified. Thus it can be concluded that a more 
sophisticated evaluation of the behaviour of the off-shell ai(1260) width should not imply a 
major qualitative and quantitative changes in our global analysis. 





a = 2 


a = 2.5 


a = 3 


Ao 

Ma 

ra,(Ml) 

xVd-o.f. 

BR(r~ 7t~^7i~tt~i't-) 


12.7 ±0.4 
1.242 ± 0.004 
0.55 ±0.02 
83.1/53 
9.11% 


11.9 ±0.3 
1.203 ±0.003 
0.48 ±0.02 
64.6/53 
9.13% 


11.9 ±0.3 
1.170 ±0.002 
0.44 ±0.01 
100.1/53 
9.21% 



Table 1: Results of three-parameter fits to ALEPH data for fixed values of a. Ma and 
r(ji(M^) are given in units of GeV. 
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Q2 (GeV2) 

Figure 3: Fit to the ALEPH data for the normahzed -K^-K^-K^Vr spectral function. 

To conclude this analysis, we take a closer look to the low region of the spectrum. 
In fact, in our approach we have assumed that corrections arising from chiral logs 

are small, hence the dominant contributions to hadronic amplitudes arise from resonance 
exchange. Close to threshold (i.e. for well below My) one is able to explicitly calculate 
the contributions of 0{p'^) chiral logs, therefore our theoretical prediction for the spectrum in 
this region can be improved, and the impact of chiral logs can be numerically evaluated. In 
order to add these contributions to the hadronic amplitude, we have made use of the calcula- 
tions performed in Ref. within standard xPT. The corresponding results are sketched in 
Appendix A. The effect on the t" — > tt+tt^tt^z/^ spectral function is shown in Fig. |5 where 
we display both the experimental values from ALEPH and the curves for the spectrum with 
and without the inclusion of O(p^) chiral logs. It can be seen that the corrections produce 
a slight enhancement of the spectral function. However, given the size of the experimental 
errors, the quality of the agreement with ALEPH data is unchanged. 

When enforcing the asymptotic conditions of QCD in Section 4, we pointed out that the 
relation in Eq. is not necessarily true when the vertices in C\^^ are considered, because 
the latter contribute to the axial form factor. If we do not take into account that constraint 
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we can perform a 5 parameter fit to the ALEPH data that gives : 

11.0 ±1.7 , 
2.55 ±0.15 , 
(1.19 ±0.02) GeV , 
(0.47 ±0.02) GeV , 

0.49 ±0.03 , (48) 

and x^/d-o.f. = 63.7/51. We see that the ratio GyjFy is consistent with the value given by 
Eq. within the error. Hence we confirm that the experimental data favour the Fy = 2Gv 
constraint. 

Finally a word of caution has to be said about the value of Aq. The only dependence 
of the amplitude on this parameter is given by the function F{Q'^,s) in Eq. ()28|) where, as 
can be seen, Aq is multiplied by the squared mass of the pion. Hence the sensibility of the 



Ao 

a 

Ma 

ra,(Ml) 

Gy 
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Figure 5: Theoretical values for the wa, wc, wd and we integrated structure functions in 
comparison with the experimental data from CLEO-II (solid) and OPAL (dashed) jH] • AH of 
them have been normalized to their respective branching ratios. 



observables in r ^ n-n-nUr on Aq is very much covered up. Tau decays into kaons would 
provide a better testing field to find out on Aq and check the results against the values from 
our fits. We consider that our result on this parameter has to be taken with care. 



5.2 Description of structure functions 

Structure functions provide a full description of the hadronic tensor T^T* in the hadron rest 
frame. There are 16 real valued structure functions in r — P1P2P3 Vj. decays (Pj is short for 
a pseudoscalar meson), most of which can be determined by studying angular correlations 
of the hadronic system. Four of them carry information on the = l"*" transitions only : 
wa, wc, wd and we (we refer the reader to Ref. [T3] for their precise definitions). Indeed, 
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for the r ttttttz/t- processes, other structure functions either vanish identically, or involve 
the pseudoscalar form factor Fp, which appears to be strongly suppressed above the very 
low-energy region due to its proportionality to the squared pion mass. 

Both CLEO-II and OPAL have measured the four structure functions quoted above for the 
— s> 7r^7r°7r°z/,- process, while concluding that other functions are consistent with zero within 
errors. Hence we can proceed to compare those experimental results with the description that 
provides our theoretical approach. In our expressions for the structure functions we input 
the values of the parameters obtained from the fit in Eq. ()47|). getting the theoretical curves 
shown in Fig. El The latter are compared with the experimental data quoted by CLEO and 
OPAL jH]. For wc, wd and we, it can be seen that we get a good agreement in the low 
region, while for increasing energy the experimental errors become too large to state any 
conclusion (moreover, there seems to be a slight disagreement between both experiments at 
some points). 

On the other hand, in the case of the integrated structure function wa, the quoted experi- 
mental errors are smaller, and the theoretical curve seems to lie somewhat below the data for 

^ 1.5 GeV^. However, it happens that wa contains essentially the same information about 
the hadronic amplitude as the spectral function dT/dQ^. This becomes clear by looking at 
Eq. (fT^ if the scalar structure function Wsa is put to zero (remember that it should be 
suppressed by a factor 0{ml/Q'^)). Taking into account that wa is given by 

wa{Q^) = JdsdtWA{Q\s,t) , (49) 

where Wa is the structure function previously introduced in Eq. (pUj) . one simply has 

dT Gp iKidp 



dQ^ 384(27r)5M, 





wa{Q') . (50) 



In this way one can compare the measurements of wa quoted by CLEO-II and OPAL 
with the data obtained by ALEPH for the spectral function, conveniently translated into 
Wa- This is represented in Fig. IHl where it can be seen that some of the data from the 
different experiments do not agree with each other within errors. Notice that, due to phase 
space suppression, the factor of proportionality between wa{Q^) and dV/dQ"^ in Eq. (jSU)) 
goes to zero for — >■ M^, therefore the error bars in the ALEPH points become enhanced 
toward the end of the spectrum. Notwithstanding, up to ^ 2.5 GeV^, it is seen that 
ALEPH errors are still smaller than those corresponding to the values quoted by CLEO-II 
and OPAL. On this basis, we have chosen here to take the data obtained by ALEPH to 
fit the unknown theoretical parameters in the hadronic amplitude. Finally, notice that a 
non vanishing contribution of Wsa (which is a positive quantity) cannot help to solve the 
experimental discrepancies, as it would go in the wrong direction. 

In the analysis of data carried out by the CLEO Collaboration [23 onto their 
7r~TT'^TT^i'r results it was concluded that the data was showing large contributions from inter- 
mediate states involving the isoscalar mesons /o(600), /o(1370) and /2(1270). Their analysis 
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T 




Q2 (GeV2) 

Figure 6: Comparison between the experimental data for wa-, from Tr^vr^vr^z/^, quoted 

by CLEO-II and OPAL and the values arising from ALEPH measurements of r~ 
71^71^71^ Uj. spectral functions jS]. The solid line is obtained from our four-parameter fit of 
the spectrum (see text). 

was done in a modelization of the axial-vector form factors that included Breit-Wigner 
functions in a Kiihn and Santamaria inspired model. Our results in the Effective Theory 
framework show that, within the present experimental errors, there is no evidence of relevant 
contributions in r — ^ TiTiTiUr decays beyond those of the p(770) and ai(1260) resonances. 

We also find interesting to study the integrated structure functions wc, wd and we, 
normalized to wa- This quantities should be less sensitive to the parameterization of the 
ai(1260) width, and do not depend on the global normalization used by each experiment 
(which could be the origin of the partial discrepancies in the values for wa)- The correspond- 
ing comparison between the theoretical and experimental results are plotted in Fig. [7[ where 
once again it is seen that our predictions are in good agreement with the experimental data 
within their present errors. 
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Figure 7: Experimental data for wc, wd and we integrated structure functions, normalized 
to wa- Solid (dashed) bars correspond to the values quoted by CLEO-II (OPAL) [H]. 
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6 Conclusions 



The hadronization of QCD currents yields relevant information on non-perturbative fea- 
tures of low-energy Quantum Chromodynamics. At -C My, xPT has provided model- 
independent knowledge on the corresponding form factors, but when resonance degrees of 
freedom become active, at ~ My, they have to be properly included into the Effective 
Theory and xPT is no longer the appropriate scheme to work with. While thorough parame- 
terizations of the hadron matrix amplitudes involving resonances have been used in the past, 
most of them implement Breit-Wigner functions in order to provide the dynamical descrip- 
tion of the resonances and, as far as we know, there is no deductive connexion between those 
functions and QCD. Hadronic r decays take place in the bulk of the energy resonance region 
and, moreover, several experiments in the last years have payed attention to them, supplying 
a good quality and quantity of data. Therefore they constitute an interesting subject of 
research to which we have turned over. 

In this article we have proposed an Effective Theory approach in order to evaluate the 
dominant form factors in r ttttttz/t- decays, and our results have been compared with the 
experimental figures. To proceed to the construction of the relevant Effective Theory we have 
relied in the R^T proposed in Ref. which has been improved here by implementing the 
necessary vector-axial-vector-Goldstone interactions ^^^^ in Eqs. fllll I12|) . The basic prin- 
ciples underlying this procedure are the all-important chiral symmetry of massless QCD, that 
drives the interaction of Goldstone bosons at low-energy, and the unitary flavour symmetry 
of the octets of resonances. Once specified our interaction Lagrangian we have proceeded to 
evaluate the axial-vector current form factors in r ^ ttttttz/t- processes. It is clear though 
that, up to this point, the Lagrangian -Cr^t is not a proper Effective Theory of QCD for 
arbitrary values of the coupling constants. Hence we have implemented the consequences 
of the asymptotic behaviour of vector and axial-vector spectral functions of the underlying 
QCD into the corresponding form factors. This procedure has fixed up several constraints 
onto the couplings, yielding an Effective Theory more germane to QCD. 

After applying this method we are left with four unknown parameters in our results for 
the axial-vector form factors. We have used the ALEPH data on r~ — *• Tr"'"TT~Tr~i/T- to perform 
a fit to its branching ratio and spectral function, rendering values for the mass and on-shell 
width of the ai(1260) resonance. Once the parameters in the theoretical expressions of the 
form factors have been determined, we have evaluated the integrated structure functions wa, 
wc, wd and we measured by the CLEO-II and OPAL experiments 0. We have found a 
good description of data, though there appear to be some inconsistencies between the results 
provided by the different experiments. 

In summary we have provided an Effective Theory based evaluation of the axial-vector 
form factors in r ^ ttttttz/,- decays that gives an appropriate account of the main features of 
the experimental data. The procedure does not rely in any modelization of the form factors, 
as it has been done in the past, but in a field theory construction that embodies the relevant 
features of QCD in the resonance energy region, showing that this is a compelling framework 
to work with. 
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A The 0{p^) xPT result for r' n+n-Jv-Ur 

In the very low region (typically E -C My), our results for the hadronic amplitudes in 
r — h'K'kVt can be compared with those obtained in standard xPT. Here we will rely on the 
analysis in Ref. where these calculations have been performed in detail up to 0{p'^). 

The O(p^) hadronic matrix element involves renormalized coupling constants ^[(/x), 
which have to be determined experimentally at a given renormalization scale. It is seen that 
these couplings can be related to finite and scale-independent quantities li j2j according to : 

^''-(/i)-ln^, (A.l) 



327rV fi 

where the coefficients 7^ arise from the corresponding renormalization group equations. In the 
resonance chiral effective theory, the constants ii are assumed to be saturated by resonance 
exchange, which at low energies induces the local xPT Lagrangian of 0{p^) for the light 
pseudoscalar mesons ^|. In this way, after considering the relations ()34|) . imposed by QCD, 
one gets 

i. = _487r2— -In 



M2 

247r2— r - In 



M2 M2 



, ml 
— m 



F 



967r2^-ln-^ . (A.2) 



Now, from the expressions quoted in Ref. 0], the = 1"*" piece of the hadronic amplitude 
for the decay — > n^n^n^h'r is found to be 



rpxPT 



2V2^^ 3s 



1+ 3F '^^] 2M2 ' 327r2F2 



2 

-s + 3t-M-3m^-(s + ml) F(s/m^) 
+ (3 1 - 2 ) F{t/ml) - {u - 2 ml) F{u/ml) 
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+ {s — 3t + u — ml.) In 




2\/2 



( 



(A.3) 



3F 



where 



F{x) = a In 



l + a 



(A.4) 



and the kinematical invariants s and t are defined as in Section 3. The first two terms in the 
curly brackets are those exphcitly quoted in Eq. (jl^j) . 
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